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1.  Overview 


This  grant  has  enabled  us  to  purchase  three  major  pieces  of  equipment:  1)  a  nuclear 
magnetic  resonance  spectrometer  system  for  studies  of  the  basic  physics  of  hyperpolarized 
129Xe  and  3He  gases;  2)  a  9.4  T  superconducting  magnet  with  a  3  inch  room  temperature 
bore;  3)  a  Verdi  diode-pumped  Nd:YAG  laser  to  replace  the  very  expensive  argon  ion  laser 
we  have  traditionally  used  for  pumping  our  Ti:sapphire  tunable  laser.  This  new  equipment 
has  greatly  improved  the  research  productivity  of  our  laboratory. 

This  equipment  supplements  our  AFOSR-sponsored  research  on  the  fundamental 
physics  and  applications,  of  spin  polarized  atoms.  One  of  the  most  unexpected  appli¬ 


cations  of  Me,  b^t£^en^in^^cal ^  miaging  of  ( hurpan  lungs ,  with  laser-polarized 
3He  and  129Xe  gas.  Basic  research,  supported  by  AFOSR  in  our  group,  eventually  led  to 


ways  to  produce  liter-atmospheres  of  3He  and  129Xe  gases  with  spin  polarizations  of  many 
tens  of  percent,  instead  of  the  parts  per  million  that  is  normal  for  magnetic  resonance 
imaging.  When  inhaled  these  gases  give  spectacular  images  of  human  lungs,  sinuses  and 
other  air  -passages.*  ‘Sphere  is'a  vigorous  'Clinical  “center  for  lung  imaging  with  hyp^rpo- 
larized  gases  at  the  University  of  Virginia  in  Charlottesville. .  That  work  is  sponsored  by 
the  international  health  care  company,  Nycomed  Amersham,  which  purchased  Magnetic 
Imaging  Technologies,  Inc.,  a  start  up  company  formed  a  few  years  ago  to  commercialize 
this  AFOSR  technology.  One  of  the  most  interesting  gas-imaging  activities  underway  now 
is  a  long-term  study  of  the  lungs  of  New  York  City  firemen  who  survived  the  collapse  of 
the  World  Trade  Center  towers  on  September  11.  2001. 

The  work  supported  by  this  AFOSR  grant  has  had  and  will  continue  to  have  appli¬ 
cations  to  many  technologies  of  the  Unites  States  Air  Force.  -  example,  tire  Rb  clocks 


used  in  the  GPS  satellite  system  are  based  on  optical  pumped  gas  cells.  The  basic  physical 
processes  operating  in  these  clocks  is  closely  related  to  that  of  the  cells  used  to  polarize  3He 
and  129Xe  gas  for  medical  imaging.  During  the  research  period  covered  by  this  report,  we 
accidentally  discovered  that  there  are  certain  novel  magnetic  resonance  transitions,  suit¬ 
able  for  use  in  atomic  clocks,  that  are  almost  unaffected  by  the  spin-exchange  broadening 
that  affects  the  traditional  0-0  clock  transition. 
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2.  Disposition  of  Equipment 

30-310  MHz  NMR  Spectrometer:  Using  the  funds  from  the  grant,  we  have  assembled  a 
rugged,  state-of-the  art  NMR  spectrometer,  continuously  covering  the  range  of  frequencies 
from  about  30  MHz  to  310  MHz.  Combined  with  the  9.4  Tesla  magnet  also  purchased  with 
funds  from  this  grant,  the  spectrometer  can  be  used  to  investigate  the  magnetic  resonance 
spectra  of  nuclei  ranging  from  14N  and  33S,  with  very  low  magnetic  resonance  frequencies, 
to  3He,  with  a  very  high  frequency.  Most  importantly,  the  spectrometer  can  be  used 
to  investigate  both  magnetic  isotopes  of  xenon, 129Xe  and  131Xe.  Unique  features  of  our 
apparatus  include  a  high-power  amplifier  (model  3304C  Herley-AMT),  with  400  Watts 
of  pulsed  power  over  the  entire  frequency  range  and  a  9-digit,  ultra  stable,  computer- 
controlled  radio-frequency  source  (model  310  PTS),  with  a  16-bit,  high-resoultion  data 
acquisition  system.  The  spectrometer  is  sensitive  enough  to  detect  deuterium  nuclei  in 
ordinary  water.  More  importantly,  it  can  also  detect  the  nuclear  magnetic  resonance 
signals  of  129Xe  or  3He,  at  atmospheric  pressure  and  at  room  temperature.  These  are  the 
only  gases  of  practical  interest  for  hyperpolarization. 

9.4  T  Oxford  Instruments  Magnet:  This  magnet  is  an  outstanding  addition  to  our 
experimental  capabilities.  It  has  superb  homogeneity  (10~7  over  a  one  inch  sphere)  and 
stability  (10~7  per  day).  It  was  installed  in  record  time,  only  one  month  after  purchase.  It 
is  located  in  a  special  part  of  our  laboratory  from  which  we  exclude  magetic  tools  or  other 
equipment  that  could  be  sucked  into  the  magnet  and  damage  it.  This  new  magnet  was 
essential  for  performing  the  work  described  in  the  attached  paper,  the  first  experiment  of 
its  kind,  and  an  impprf ant  contribution  to  bur  understanding  of  spin  physics  in  gaseous 
media. 

Verdi  Pump  Laser:  This  solid-state  Nd:Yag  laser,  with  diode-laser  pumping,  has  been 
used  to  replace  an  argon-ion  laser  that  is  the  source  of  optical  pumping  light  for  our 
.  TkSapphire  laser.  The  argon-ion. laser  that  it  replaced  required  high-pressure,  deionized 
Cooling  water  from  the  university  chilled  water  plant.  Also  the  very  expensive  plasma'tubes 
needed  to  be  replaced  every  few  years.  The  new  Verdi  laser  needs  no  external  cooling  water; 
since  it  is  much  more  energy  efficient  than  the  argon-ion  laser.  A  closed  cycle  cooling  loop 
with  air  cooling,  is  all  that  Is.  required.  We  have  found  the  Verdi  laser  very  convenient  to 
use,  and  we  expect  to  save  a  great  deal  of  time  and  money  with  it  in  the  future.  The  Verdi 
laser  was  used  to  obtain  the  data  described  in  the  refereed  publication  and  in  the  attached 
preprint. 


3.  Publications 

Although  the  DURIP  equipment  was  only  installed  last  year,  it  has  already  made 
possible  new  areas  of  research  that  we  could  not  have  done  previously.  One  refereed  paper 
using  this  equipment  has  already  been  published,  and  we  attach  a  preprint  of  one  that  has 
just  been  submitted  for  review. 

D.  K.  Walter,  W.  M.  Griffith  and  W.  Happer,  “Magnetic  Slowing  Down  of  Spin 
Relaxation  due  to  Binary  Collisions  of  Alkali-Metal  Atoms  with  Buffer-Gas  Atoms,”  Phys. 
Rev.  Letters  88  093004  (2002). 
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4.  Personnel 


The  following  personnel  have  used  the  equipment  purchased  with  Grant  AFOSR 
F49620-01-1-0254: 

Dr.  William  Happer,  Principal  Investigator. 

Dr.  Nicholas  Kuzma,  Dicke  Fellow  and  Research  Associate. 

Dr.  Warren  Griffith,  Research  Associate 

Dr.  Dan  Walter,  Research  Associate,  PhD  with  AFOSR  support  in  2001. 

Mr.  Yuan  Yu  Jau,  Graduate  Student  and  PhD  candidate. 

Ms.  Amber  Post,  Graduate  Student  and  PhD  candidate. 
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High-field  measurement  of  the  129Xe-Rb  spin-exchange  rate  due  to  binary  collisions 

Yuan- Yu  Jau,  Nicholas  N.  Kuzma,  and  William  Happer 
Department  of  Physics,  Princeton  University ,  Princeton,  New  Jersey  08544 

(Dated:  July  12,  2002) 

We  have  measured  the  binary  spin-exchange  rate  coefficient  for  collisions  between  Rb  and  129 Xe 
atoms  at  a  magnetic  field  B  =  9.4  T  and  in  the  temperature  range  160  <  T  <  200°C.  The  longitu¬ 
dinal  spin  relaxation  rate,  1/Ti,  of  129Xe  is  dominated  by  129Xe-Rb  binary  collisions.  We  measured 
1/Ti  with  the  methods  of  conventional  nuclear  magnetic  resonance,  with  no  spin  enhancement  by 
optical  pumping.  At  the  high  magnetic  field  and  large  gas  pressures  used  in  this  work,  the  con¬ 
tribution  of  129Xe-Rb  molecules  to  the  spin-exchange  rate  was  less  than  1%  of  the  binary-collision 
rate.  The  Rb  atomic  number  density  [Rb]  was  measured  directly  from  the  large  Faraday  rotation 
of  the  polarization  of  near-resonant  light  passing  through  the  sample  cells.  The  rate  coefficient  is 
k  =  d(l/Ti)/d[Rb]  =  (1.75  ±  0.12)  x  10"16  cm3/s. 


PACS  numbers:  32.80.Bx,  33.25.+k,  34.30. +h 

I.  INTRODUCTION 

Spin-exchange  optical  pumping  [1,  2]  is  becoming  a 
widely-used  method  to  produce  large  quantities  of  hyper- 
polarized  3  He  [3,  4]  and  129Xe  [5].  These  hyperpolarized 
gases,  where  the  nuclear  spin  polarization  is  much  larger 
than  its  thermal-equilibrium  value,  have  found  applica¬ 
tions  in  magnetic  resonance  imaging  [6-8],  neutron  spin 
structure  studies  [9],  surface  analysis  [10,  11],  and  experi¬ 
mental  tests  of  fundamental  symmetries  [12, 13].  Many  of 
these  applications  would  benefit  from  an  increase  in  spin 
polarization  of  the  gas.  Typical  polarizations  of  129Xe. 
axe  still  a  factor  of  three  lower  than  those  of  3He. . 

Currently,  Rb  is  most  commonly  used  to  transfer  an¬ 
gular  momentum  from  the  optical-pumping  laser  light  to 
the  nuclei  of  noble  gas*  atoms.  Compared  to  the  spin- 
exchange  processes  in  the  3He-Rb  system  [14],  where  the 
exchange  occurs  almost  exclusively  in  binary  collisions, 
spin  exchange  in  129Xe-Rb  [15]  has  contributions  from 
both  binary  collisions  and  from  the  three-body  forma¬ 
tion  and  collisional  breakup  of  129Xe-Rb  van  der  Waals 
molecules.  A  better  understanding  of  the  relative  impor¬ 
tance  of  binary  collisions  and  van  der  Waals  molecules 
for  polarizing  129Xe  will  facilitate  the  design  of  systems 
that  produce  more  highly  polarized  129Xe.  In  this  pa¬ 
per  we  describe  a  new,  direct  way  to  measure  the  binary 
spin-exchange  rate  that  is  free  of  many  of  the  difficulties 
that  beset  earlier  measurements. 

In  their  pioneering  studies  of  the  129Xe-Rb  system. 
Volk  et  ai  [16]  reported  a  rate  coefficient  k,  =  45  x 
10“16cm3/s.  High-permeability  shields  were  used  to  re¬ 
duce  the  magnetic  field  to  less  than  10-9  T,  and  the 
temperature  was  T  ^  60°  C.  The  binary  rate  coefficient 
reported  by  Volk  et  al.  is  more  than  a  factor  of  ten  larger 
than  that  reported  in  most  subsequent  work.  We  men¬ 
tion  a  few  possible  reasons  for  this  large  discrepancy.  The 
experiments  of  Volk  et  al.  were  done  with  glass  cells,  con¬ 
taining  excess  Rb  metal  and  sealed  off  at  room  temper¬ 
ature  with  half  a  Torr  of  Xe  and  N2  at  partial  pressures 
ranging  between  10  and  100  Torr.  At  such  low  fields  and 
gas  pressures,  most  of  the  spin  relaxation  of  129Xe  is  due 


to  XeRb  molecules.  It  is  difficult  to  extract  the  relatively 
small  contribution  due  to  binary  collisions.  The  difficul¬ 
ties  were  compounded  by  the  oversimplified  model  used 
to  analyze  the  data.  For  example,  not  accounting  for  the 
slowing-down  of  the  relaxation  due  to  the  angular  mo¬ 
mentum  stored  in  the  nuclear  spins  of  Rb  atoms  would 
cause  one  to  substantially  underestimate  the  contribution 
of  XeRb  molecules  to  the  relaxation  rates. 

In  subsequent  experiments,  Cates  et  al  [17]  reported  a 
binary  rate  coefficient  of  k  —  (3.7  ±  0.7)  x  10*"16  cm3/s; 
about  12  times  smaller  than  that  of  Volk  'et  al.  The  mea-'* 
surements  were  done  at  temperatures  of  about  100°  C. 
oThe  magnetic  field  was  1.1  x  10~5  T,  too  small  to  cause 
appreciable  suppression  of  the  relaxation  due  to  129Xe- 
Rb  molecules.  The  model  used  by  Cates  et  al  to  analyze 
their  data  took  into  account  the  nuclear  spins  of  the  Rb 
atoms.  To  eliminate  the  effects  of  129Xe-Rb  molecules, 
they  made  measurements  in  a  series  of  four  cells,  each 
containing  (at  room  temperature)  about  50  Torr  of.nitror 
gen  and  with  xenon  pressures  ranging  from  245  to  1817 
Torr..  An  extrapolation  procedure  was  used  to  infer  the 
rate  coefficient  at  the  infinite  pressure  limit  ,  where  molec-  • 
ular  contributions  would  be  . negligible.  However,  recent 
findings  of  Kadlecek- et  al  [18],  confirmed  by  work  in'  our- 
laboratory,  show  that  the  models  of  relaxation  due  to  the 
formation  and  breakup  of  loosely  bound  molecules  ha.ve 
serious,  and  still  poorly  understood  flaws,  at  gas  densities 
on  the  order  of  one  amagat  or  greater.  This  may  have 
introduced  some  problems  in  the  pressure  extrapolation 
procedure  of  Cates  et  al  [17]. 

The  first  high-field  measurements  of  rate  coefficients 
were  reported  by  Augustine  et  al  [19,  20].  At  a  field 
of  B  —  2.35 T  and  a  xenon  partial  pressure  of  450  Torr, 
they  reported  k  =  2.8  x  10“ 1G  cm3/s.  At  the  same  field 
but  at  a  xenon  partial  pressure  of  1500  Torr,  they  re¬ 
ported  ac  =  7.4  x  10“ 16  cm3/s.  It  is  not  clear  what 
might  have  caused  this  puzzling  dependence  on  xenon 
pressure.  One  would  expect  a  decreasing  contribution  of 
XeRb  molecules  to  the  relaxation  rates  with  increasing 
pressure;  the  increasing  three-body  formation  rate  does 
not  compensate  for  the  more  rapidly  decreasing  spin  flip 
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probabilities  that'  result  from  the  shorter  molecular  life- 
times.  At  pressures  comparable  to  those  used  by  Augus¬ 
tine  et  al ,  but  at  much  smaller  magnetic  fields,  the  rate 
coefficients  measured  by  Cates  et  al  [17]  decreased  with 
pressure. 

In  all  three  of  these  early  measurements,  the  number 
density  of  Rb  atoms  [Rb]  was  not  measured  directly  but 
was  inferred  from  the  measured  temperature  T  of  the 
sample  cells.  Reference  was  made  to  empirically  deter¬ 
mined  formulae  for  the  dependence  of  [Rb]  on  T.  The 
Killian  [21]  formula  used  by  Cates  et  al  [17]  gives  val¬ 
ues  for  [Rb]  that  are  approximately  twice  as  large  as 
the  values  given  by  the  Smithells  formula  [22],  used  by 
Volk  et  al  [16]  and  by  Augustine  et  al  [19,  20].  Ex¬ 
perience  has  shown  that  different  cells  containing  Rb 
metal  at  the  same  temperature  can  have  Rb  vapor  den¬ 
sities  that  would  differ  by  factors  of  two  or  three  from 
the  formulas  of  Killian  or  Smithells.  The  early  experi¬ 
ments  [16,  17,  19,  20]  were  done  at  relatively  low  [Rb] 
values,  typically  [Rb]  <  1011  cm-3.  Contributions  of  the 
cell  walls  [23],  magnetic-field  gradients,  and  other  mecha¬ 
nisms  [24,  25]  to  the  relaxation  rates  of  129Xe  could  have 
affected  these  earlier  measurements. 


II.  EXPERIMENT 

Here  we  ,  report  the  first  direct  measurement  of  the 
129Xe-Rb  spin-exchange  rate  due  to  binary  collisions.  We 
chose  high  .gas  pressures  (>  600  Torr)  and  a  high  ipag- 
:  infetic  feeld  B  9%  T  Sb  th&t  the  dontfiblif ion"  of  RbXe 
molecules  to  the  spin-exchange  rate  would  be  negligible 
compared  to  the  contribution  of  binary  collisions.  We 
used  nuclear  magnetic  resonance  (NMR)  to  probe  the 
thermal-polarization  signal  of  129Xe  nuclei  without  opti¬ 
cal  pumping,  and  we  used  the  Faraday  rotation  method 
[26]  to 'directly  determine  [Rb].' All  measurements  were 
done  at  high  temperatures,  between  160°  and  200°  C, 
and  consequently  at  much  higher  Rb  densities  than  in 
the  previous  studies  [16,  17,  19,  20]. 

The  return  of  the  129 Xe  nuclear  spin  polarization  (if*) . 
to  its  thermal  equilibrium  value  (KZ)T  =  pLxB/2kT 
(here  k  is  the  nuclear  magnetic  moment  of  129Xe  and  k 
is  the  Boltzmann  constant)  after  saturation  (destruction) 
by  a  train  of  radio-frequency  (rf)  pulses  is  described  by 

(Kz)  =  (Kz)t( l-e-t/T>)  ,  (1) 

where  we  can  write  the  longitudinal  relaxation  rate  1/Ti 
of  the  129Xe  nuclei  as 

i-  =  [Rb]«+r.  (2) 

In  this  expression,  k  (also  denoted  as  (aex^)  in  the  lit¬ 
erature)  is  the  rate  coefficient  for  binary  spin-exchange 
collisions  between  129Xe  and  Rb  atoms,  and  T'  is  the 
spin-relaxation  rate  due  to  other  causes,  such  as  contri¬ 
butions  from  129Xe-Rb  van  der  Waals  molecules  [1],  wall 


interactions  [23],  the  formation  of  Xe2  molecules  [24],  or 
binary  Xe-Xe  collisions  [25].  We  will  present  both  theo¬ 
retical  (App.  A)  and  experimental  evidence  that  the  val¬ 
ues  of  1/Ti  measured  in  this  work  are  dominated  by  the 
binary  collision  contribution,  [Rb]  rc. 

In  our  work,  we  measure  the  longitudinal  129Xe  relax¬ 
ation  rate  1/Ti  at  different  cell  temperatures  T.  The 
main  effect  of  changing  the  cell  temperature  is  the  large, 
exponential  change  in  Rb  density  [Rb]  a  exp(-T4/T) 
[the  Arrhenius  temperature  «  9500  K  »  T,  so  that 
[Rb]  oc  exp(TAAT/T2)  «  exp(AT/20K)  ].  The  temper¬ 
ature  dependence  of  the  other  relaxation  rates  T'  is  ex¬ 
pected  to  be  so  small  that  we  take  dV /dT  =  0  over  the 
temperature  range  of  our  measurements.  Our  experimen¬ 
tal  data  are  consistent  with  neglecting  any  temperature 
variation  of  T'.  As  we  mentioned  before,  the  cell  temper¬ 
ature  T  does  not  uniquely  define  [Rb],  so  at  each  value 
of  T  we  independently  determined  [Rb]  by  measuring  the 
Faraday  rotation  of  the  plane  of  polarization  of  a  prob¬ 
ing  laser  beam  passing  through  the  cell.  Our  modeling 
calculations  indicate  that  there  should  be  little  change 
in  the  binary  exchange  rate  coefficient  within  the  range 
160  <  T  <  200°  C. 

The  mean  binary  spin-exchange  rate  coefficient,  k,  is 
therefore  the  slope  of  the  experimentally  determined  plot 
of  1/Ti  versus  [Rb].  The  two  quantities,  1/Ti  and  [Rb], 
must,  be  determined  with  good  experimental  accuracy 
over,  a  range  of  cell  temperatures.  . 

A.  Sample  Preparation 

In  this  work,  we  made  four  spherical  .aluminp-silicate 
*  glass  cells,  2.5  cm  in  diameter.  These  were  attached  to 
a  glass  tube,  connected  at  one  end  to  a  gas-handling 
manifold  with  pressure  .mqnfrormg  equipment  and.  with 
cylinders  of  ultra-pure  nitrogen  and  isotopically  enriched 
xenon  (71%  129Xe).  The  manifold  could  be  evacuated 
with  a  turbo- vacuum  pump.  Attached  to  the  other  end 
of  the  glass  tube  was  a.  glass  retort  with  a- mixture  of 
RbCl  salt  and  chips  of  Gai  metal.  The  cells  were  baked 
overnight  at  400°  C  under  a  vacuum  of  10“ 7  Torr.  Then  a 
gas  flame  was  used  to  heat  the  retort  and  release  volatile 
Rb  metal.  A  few  milligrams  of  rubidium  metal  were  dis¬ 
tilled  into  each  of  the  cells  with  a  gas  torch.  A  cell  was 
then  immersed  in  liquid  nitrogen  to  condense  all  of  the 
xenon  gas  of  known  room-temperature  pressure  and  vol¬ 
ume.  Nitrogen  gas  at  less  than  atmospheric  pressure  (so 
it  could  not  condense  in  the  cold  cell)  was  added  to  the 
manifold  and  the  cell  stem  was  sealed  off  while  the  body 
of  the  cell  remained  immersed  in  liquid  nitrogen.  With 
the  ideal  gas  law  and  the  known  volume  of  the  cell,  it  was 
possible  to  determine  the  number  densities  [Xe]  and  [N2], 
which  we  quote  in  units  of  amagats  (1  amg  ~  2.69  x  10 19 
cm-3).  We  prepared  cells  with  four  different  densities  of 
nitrogen  and  xenon.  After  being  sealed-off,  each  cell  was 
put  into  an  oven  at  200°  C  for  at  least  three  days  in  order 
to  make  the  wall  relaxation  of  129Xe  reach  a  stable  value 
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FIG.  1:  Schematic  diagram  of  the  experimental  apparatus 
built  to  measure  1  /7\  and  [Rb]  for  four  sealed-off  sample  cells. 
More  details  are  given  in  the  text.  The  rf  frequency  of  the 
NMR  spectrometer  was  offset  from  the  resonant  129Xe  fre¬ 
quency  by  ~  1  kHz  in  order  to  avoid  the  low  frequency  noise 
background.  The  wavelength  scanning  range  was  from  765  to 
795  nm.  The  oven  was  ohmically  heated  up  to  200°  C. 


in  the  presence  of  rubidium  through  a  “curing  process” 

in].  . .  :  .  .. 

B.  NMR  Measurement  of  129Xe  Relaxation  Times 

Our  experimental  setup  is  shown  in  Fig,  1.  The  cell 
and  a  surrounding  NMR  coil  were  mounted  inside  a 
heat-insulated,  cylindrical  oven,  suspended  in  the  89- 
mm  room-temperature  bore  of  a  9.4-T  Oxford  Instru¬ 
ments  superconducting  magnet.  The  inner  volume  of  the 
oven  was  heated  ohmically  using  weakly-magnetic  high- 
resistance  wire,  counter-wrapped  8  cm  below  the  sample 
cell  to  avoid  distorting  the  3  xlO-7  homogeneity  of  the 
magnetic  field.  The  temperature  was  monitored  and  con¬ 
trolled  with  two  resistive  platinum  sensors  (RTD)  and  an 
Omega  controller.  For  the  T\  measurement,  a  previously 
described  NMR  spectrometer  [27]  was  used  at  a  radio  fre¬ 
quency  set  to  110.4800  MHz.  The  pulse  sequence  started 
with  a  saturating  pulse  train  (100  pulses  500  //s  each)  that 
was  sent  to  destroy  the  thermal  polarization  of  129Xe. 
After  a  variable  recovery  time  rr,  a  single  90°  pulse  was 
used  to  excite  a  free  induction  decay  (FID)  signal.  By 
varying  rr,  we  acquired  many  saturation-recovery  curves 
like  the  one  shown  in  Fig.  2. 


C.  Faraday  Rotation  Measurement  of  Rb  Density 

The  rubidium  atomic  number  density  [Rb]  was  directly 
determined  by  a  Faraday  rotation  measurement.  A  Ti- 
sapphire  ring  laser  “Coherent  899-29”,  pumped  by  a  Nd- 


FIG.  2:  A  representative  saturation-recovery  measurement 
of  the  longitudinal  relaxation  time  T\  of  129 Xe  nuclear  spins 
at  a  given  temperature.  Details  of  the  NMR  procedure  are 
described  in  the  text.  Each  data  point  is  the  average  of  many 
FID  signals  with  the  same  recovery  time. 


YAG  “Verdi”  laser,  was  used  to  generate  the  probing 
beam  and  to  scan  the  wavelength  of  light.  We  used  a 
linear  polarizer  (LP)  to  ensure  that  the  light  entering 
the  cell  was  linearly  polarized. 

We  can  describe  a  light  beam  in  a  pure  polarization 
state  with  two  parameters:  (i)  the  orientation  angle  9 
of  the  semimajor  axis  of  the  elliptically  polarized  light 
relative  to  the  a;  axis  in  a  reference  plane  normal  to  the 
beam  direction,  and  (ii)  the  mean  photon  spin  s,  a  real 
number  ranging  from  ^-1  to  4*1,  which  describes  the  rel¬ 
ative  circular  polarization  of  the  light.  The  ratio  a  of  the 
minor  to  the  major  axis  of  the  ellipse  is  related  to  the 
mean  spin  s  by  a  =  (1  —  y/l  —  s2)/s. 

In  our  experiment,  the  semimajor  axis  of.  the  ellip¬ 
tically  polarized  light  entering  the  cell  was  oriented  at 
an  angle  9q,  which  we  could  vary  for.  experimental  coin 
venience.  The  mean  spin  so  of  the  light  entering  the 
cell  was  close  to  zero,  since  the  probe  light  had  nearly 
perfect  linear  polarization.  Stray  birefringence  in  lenses 
and  windows  and  reflections  from  two  turning  mirrors, 
not  shown  in  Fig'.  1,  could  change  both  0  and  s,  but  we 
verified  experimentally  that  these  were  negligible  effects. 
The  orientation  angle  of  the  light  emerging  from  the  cell 
is  0  —  Qq  4-  A0,  where  the  Faraday  rotation  angle  A 0 
can  be  several  radians  because  of  large  [Rb]  and  B.  In 
principle,  the  circular  dichroism  of  the  vapor  can  induce 
a  corresponding  change,  As,  in  the  mean  photon  spin, 
but  unless  the  laser  wavelength  is  very  close  to  the  cen¬ 
ter  of  the  absorption  lines,  As  will  be  negligible  and  the 
light  emerging  from  the  vapor  will  remain  nearly  linearly 
polarized  with  s  ~  0. 

To  measure  the  Faraday  rotation  angle  A9 ,  we  set  the 
stress  axis  of  a  photo-elastic  modulator  (PEM)  at  45° 
to  the  axes  of  a  linearly-polarizing  beam-splitter  cube 
(BS-LP)  [28].  The  angle  rj)  of  the  relative  phase  re¬ 
tardation  between  the  fast  and  slow  axes  of  the  PEM 
was  varied  sinusoidally  with  time  t  at  the  frequency 
o;p/27r  =  83.4  kHz.  That  is,  ip  =  ^sincupt,  where  the 
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peak  phase  shift  tpp  is  proportional  to  the  driving  volt¬ 
age  of  the  PEM  and  inversely  proportional  to  the  laser 
wavelength.  After  passing  through  the  BS-LP,  the  probe 
beam  was  detected  with  a  photodiode,  which  produced 
the  output  voltage, 


where  the  constant  G  depends  on  the  different  LIA  gains 
in  the  system.  From  the  well-known  theory  of  Faraday 
rotation  (e.  g.  see  [26]),  the  rotation  angle  for  the  light 
of  wavelength  A  =  2iccjo)\  close  to  the  strong  and  D2 
resonance  lines  is 


Vpd  =  glB(  l  +  \/l-s2  sin  26  cos  ±  s  sin  -0),  (3) 

where  g  is  the  “gain”  of  the  detection  system.  The  laser 
power  Iq  =  1#(£)  was  mechanically  chopped  at  the  fre¬ 
quency  ft/27r  =  153  Hz. 

The  laser  wavelength  A  was  scanned  with  time,  and 
this  caused  large  variations  in  the  angle  0,  and  also  some 
variation  of  Ib  because  of  the  wavelength  dependence  of 
the  laser  output  power.  There  were  also  noise  fluctua¬ 
tions  of  /£. 

We  can  decompose  the  signal  (3)  into  different  fre¬ 
quency  harmonics  in  the  familiar  way  for  the  frequency 

modulation,  Vpd  =  Vo  +  K,p  smupt+ V2Wp  cos  2 upt  H - , 

where  amplitudes 

Vo  =  gls  *  (1  +  \/l  -  52  J0(^p)  sin  26)  ,  (4) 

=  2 gls  *  sJi(i/)p)  ,  (5) 

V2u,p  =  2 gls  •  \/l  -  s2J2(V'p)  sin  20  (6) 

are  chopped  at  the  relatively  low  frequency  SI.  Here  Jn 
are  Bessel  functions  of  the  first  kind.  We  empirically 
adjusted  the  PEM  driving  voltage  to  find  the  smallest 
amplitude  i>p  =  (that  is,  *00  =  2.405)  for  which 
J0(t/>p)  =  0.  Then  Vo  =  gls  became  independent  of 
the  linear  polarization  angle  0,  and  could  be  used  to  nor¬ 
malize  the  signals  VWp  and  V^p .  A  small  correction  was 
made  to  the  data-to -account  for-the  slight  ^variation  of 
ij)p  oc  1/A  over  the  1%  range  of  the  wavelengths  used  in 
this  work. 

As  sketched  in  Fig.  1,  the  signal  V^p  was  detected 
with  two  lock-in  amplifiers  (LIA).  The  voltage  Vpd  was 
applied  to  the  first  LIA,  which  was  referenced  to  the  PEM 
drive  frequency  2up  with  a  1  ms  integration  time.  The 
output  of  the  first  LIA  was  applied  to  a  second  LIA,  ref¬ 
erenced  to  the  probe  laser  chopping  frequency  SI  with  a 
0.1  s  integration  time.  For  further  analysis,  the  output 
voltage  of  the  second  LIA  is  processed  with  one  channel 
of  a  digitizing  interface  (not  shown  in  Fig.  1)  to  obtain  a 
digitized  signal  S2ujp  oc  V^p*  This  was  stored,  along  with 
the  laser  wavelength,  in  a  computer  (PC).  The  signal  Vo 
was  detected  by  applying  Vpd  to  a  third  LIA,  referenced 
to  the  chopper  frequency  SI  and  with  a  0.1  s  integration 
time,  and  acquired  as  described  above  to  obtain  the  digi¬ 
tized  signal  So  oc  Vo*  A  digitized  record  of  the  probe  laser 
wavelength  A  was  stored  along  with  each  pair  of  values, 
So  and  S2LJp. 

For  further  analysis  of  the  digitized  data,  we  computed 
the  ratio  R  =  J?(A)  =  S^/So-  According  to  (4)  and  (6) 
the  ratio  should  be 


R  =  S2up/SQ  =  ic  sin  2(0o  +  A0),  (7) 


A0  =  l  [Rb]  2 

Jmcr 


Aj  mcr  (UJX  ~  COjmo) 
[(wA  -  UjmaY  +  72/4] 


where  /  is  the  optical  path  length  of  the  laser  beam  inside 
the  cell  and  the  the  quantum  numbers  Jmcr  are  discussed 
below.  The  coefficients  Ajma  are 


Ajrna  — 


3tt u>\rec  fj  a  C('™ 


m+CT 

m;  1,  a 


Pm 


(2  J  +  1)  U)j 


(9) 


where  re  =  e2/mec2  =  2.82  x  10”13cm  is  the  classi¬ 
cal  electron  radius,  e  and  me  are  the  electron  charge 
and  mass  in  cgs  units,  c  is  the  speed  of  light,  fj  are 
the  oscillator  strengths  of  absorption  lines  to  the  52P.; 
excited  states.  To  a  good  approximation  fi/2  =  1/3 
and  /3/2  =  2/3.  The  resonant  frequencies  of  the  two 
D  lines  are  ojj  =  2nc/Xj.  These  are  pressure  shifted 
slightly  from  the  free-atom  values  of  A1/2  =  7947  A 
and  A3/2  =  7800  A.  We  assume  pressure-broadened 
Lorentzian  lines  with  the  same  width  at  half  maximum 
of  7,  as  indicated  in  Eq.  (8).  The  quantity  in  square 
brackets  is  a  Clebsch-Gordan  coefficient  [29].  The  index 
a  =  ±1  determines  whether  the  absorbed  light  is  left  or 
right  circularly  polarized.  The  index  m  =  ±1/2  is  the 
azimuthal  electron  spin  quantum  number  of  the  ground 
state  Rb  atom,  and  pm  —  Z~l  exp{-2 mnsR/kT}  is  the 
probability  of  a  given  atom  being  in  a  state  with  the 
quantum  number  m.  Here  Z  =  2  cosh  //^B/feT;  the  cell 
temperature  is  T  (in  kelvin),  and  /13  is  the  Bohr  magne¬ 
ton. 

The  resonant  frequencies  uijma  are 


Wjma  =  Uj  -b 


+  -  2m], 


(10) 


where  the  Zeeman  splitting  of  the  excited-state  az¬ 
imuthal  levels  is  determined  by  the  g- values,  gi/2  — 
2/3  and  g$/2  =  4/3.  Of  the  eight  possible  combi¬ 
nations  for  the  summation  indeces  in  Equation  (8), 
Jmcr  =  (1  ±  1/2,  ±1/2,  ±1),  the  triples  (1/2, 1/2, 1) 
and  (1/2, —1/2, —1)  are  excluded  since  they  refer  to 
non-existent  excited-state  sublevels  with  J  =  1/2  and 
mj  =  ±3/2.  So  there  are  only  6  non-zero  terms  in  the 
sum  (8). 

From  inspection  of  Eq.  (7),  we  see  that  the 
experimentally-determined  ratio  R  on  the  left  can  be 
made  equal  to  the  theoretical  expression  on  the  right  by 
setting  the  gain  parameter  G  to  twice  the  maximum  value 
of  R,  and  then  fitting  the  unknown  number  density  [Rb] 
and  the  approximately  known  initial  orientation  angle  60 
and  the  resonance-line  wavelengths  A  j. 

With  the  relatively  large  values  of  magnetic  field 
(B  =  9.4T)  and  atomic  number  densities  ([Rb]  ^ 
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FIG.  3:  Experimentally  determined  Faraday  rotation  angles 
A 6  and  the  best  fitting  curve,  based  on  (8).  Our  apparatus 
did  not  permit  us  to  record  the  negative  rotation  angles  that 
occur  at  small  detunings.  The  large  rotation  angle  A 6  of  the 
plane  of  polarization  permits  an  unambiguous  fit  of  the  data 
to  the  theory,  and  gives  an  excellent  measurement  of  the  Rb 
number  density  [Rb]. 

1014cm-3)  used  in  our  experiments,  we  could  observe 
several  oscillations  of  the  signal  V^p ,  with  negligible  at¬ 
tenuation,  as  the  laser  was  tuned  toward  the  resonance 
line.  The  resulting  highly-structured  signal  could  be  fit 
very  precisely  (Fig.  3).  Our  digitizing  interface  only  ac¬ 
cepted  voltages  in  the  range  0  -  10  V,  thus  the  negative 
data  points  are  not  shown  in  Fig.  3.  To  obtain  the  max¬ 
imum  data  with  this  limitation,  the  initial  orientation 
angle  0O  was  chosen  to  make  the  signal  from  the  second 
LI  A  slightly  positive  in  the  limit  of  large  positive  detun¬ 
ings,  —  u>3/2  7-  Our  data,  such  as  shown  in  Fig.  3, 

are  much  easier  to  fit  than  the  previous  Faraday  rotation 
data  [26,  28],  taken  at  much  lower  values  of  [Rb]  and  jB, 
where  A 8  <£.  n/4  and  thus  useful  Faraday  rotation  sig¬ 
nals  could  be  observed  only  for  detunings  much  smaller 
than  the  D-line  splittings.  In  present  work,  we  acquired 
Faraday  rotation  data  with  good  signal-to-noise  ratios  for 
detunings  comparable  to  the  D  line  splittings,  so  it  was 
necessary  to  account  for  both  Dj  and  D2  lines  in  the  ex¬ 
pression  (8),  rather  than  only  one,  which  was  appropriate 
for  the  earlier  work  [26,  28]  with  smaller  detunings. 

D.  Measurement  of  Stray  Circular  Polarization 

In  our  experiment,  dielectric-coated  mirrors  are  used  to 
change  the  initial  horizontal  direction  of  the  laser  beam 
to  vertical  for  propagation  up  along  the  bore  of  the  super¬ 
conducting  magnet.  Above  the  magnet,  a  second  mirror 
was  used  to  change  the  beam  direction  back  to  horizontal 
to  facilitate  mounting  of  the  PEM  and  the  light-detection 
apparatus.  To  check  for  any  polarization  changes  pro¬ 
duced  by  the  mirrors,  we  replaced  the  reference  signal 
at  2lop  in  the  first  LIA  of  (Fig.  1)  with  a  reference  sig¬ 
nal  at  u}p.  With  this  arrangement,  the  output  of  the 


Rubidium  number  density  [Rb]  (1014cm'3) 


FIG.  4:  Measured  longitudinal  relaxation  rates  1/Ti  plotted 
versus  measured  Rb  densities  [Rb]  for  the  four  cells  listed  in 
Tab.  I.  The  spin-exchange  rate  coefficients  ft  are  the  slopes  of 
the  straight-line  fits  to  the  data. 


Cell  No. 

Xe  density  [amg] 

N2  density  [amg] 

ft  [10  w  cnF/s] 

T097 

4.0 

0 

1.78  ±0.10 

TO  96 

0.79 

0 

1.76  ±0.16 

T095 

0.79 

1.0 

1.60  ±0.17 

TO  94 

0.79 

3.0 

1.84  ±0.17 

TABLE  I:  Names,  gas  content,  and  experimentally  measured 
values  of  k  for  the  four  sample  cells  used  in  our  experiment. 


second  LIA  was  proportional  to  VUJp  and  thus  s,  as  one 
can  see  from  Eq.  (5).  In  this  way  we  verified  |s|  <  0.07, 
so  our  turning  mirrors  introduced  negligible  phase  shifts 
between  orthogonal  senses  of  polarization  of  the  probe 
beam.  Other  dielectric  mirrors,  not  used  in  this  experi¬ 
ment,  introduced  quite  large  phase  shifts. 


III.  RESULTS  AND  DISCUSSION 

Figure  4  is  a  summary  plot  of  the  experimentally  mea¬ 
sured  1/Ti  values  versus  the  corresponding  values  of  [Rb], 
measured  with  the  Faraday  rotation  technique,  for  the 
four  cells  listed  in  Tab.  I  .  We  fit  the  data  for  each 
cell  to  a  straight  line  to  obtain  a  slope  and  an  inter¬ 
cept.  In  accordance  with  Eq.  (2),  the  slopes  are  our  es¬ 
timates  of  ft,  while  the  intercepts  are  our  estimates  of 
the  additional  relaxation  F'.  Within  the  experimental 
error,  the  four  slopes  are  identical.  This  experimen¬ 
tally  validates  our  neglecting  any  temperature  depen¬ 
dence  of  T'  in  the  earlier  analysis.  The  values  of  V 
extrapolated  from  the  data  shown  in  Fig.  4  are  con¬ 
sistent  with  r  =  0  to  within  the  experimental  errors 
(a  direct  measurement  of  F'  at  room  temperature  gave 
T'  =  0.002±0.001  s”1).  The  mean  slope  for  the  four  data 
sets  was  ft  =  (1.75  ±  0.12)  x  10“16cm3/s. 

We  have  measured  k{B,T)  at  a  magnetic  field  B  = 
9.4  T  and  temperatures  T  from  433  to  473  K.  Rased  on 
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the  modeling  calculations  we  have  done,  both  with  clas¬ 
sical  path  calculations  and  with  the  distorted  wave  Born 
approximation,  we  expect  the  zero-field,  100°  C  (373  K) 
value  k(0,  373  K)  to  be  about  20%  larger  compared  to  the 
«(9.4  T,  453  K)  value  measured  here.  The  9.4-T  magnetic 
field  is  sufficient  to  partially  suppress  the  spin-flip  transi¬ 
tion  probabilities,  even  for  binary  collisions,  which  have 
durations  of  only  picoseconds.  We  will  report  on  these 
modeling  calculations  in  a  subsequent  paper. 
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APPENDIX  A;  HIGH-FIELD  CONTRIBUTION 
OF  VAN  DER  WAALS  MOLECULES 

For  simplicity  of  notation  we  assume  100%  isotopic 
abundance  of  129Xe.  We  can  write  the  spin  Hamiltonian 
of  a  XeRb  molecule  as 

H  =  TUoSz  +  aS-K  =  HW  +HW.  (Al) 

Here,  S  is  the  spin  operator  of  the  Rb  valence  electron, 
K  is  the  spin  operator  of  the  129Xe  nucleus,  a/ fizz  2.0  x 
108  s-1  is  a  representative  value  of  the  hyperfine  coupling 
coefficient,  averaged  over  vibrational-rotational  states  of 
the  molecule,  and  u  =  2 fiBB/h  =  1.65  x  1012  s"1  is  the 
Larmor  frequency^  theeleetroniir  a  magnetkriield  of 
9.4  T.  At  this  large  field,  we  can  ignore  the  nuclear  spin 
of  the  Rb  atom  in.  calculating  its  spin  exchange  rates 
with  129Xe.  The  zeroth-order  part  of  the  Hamiltonian 
(Al)  is  H —  hu>Sz  +  cxSxKz  and  the  spin-exchange 
perturbation  is  H ^  =  a(S+K _  +S-K+)/ 2. 

The  129XeRb  molecule  will  be  formed  in  one  of  the 
four  spin  states  | where  ms  —  ±1/2  and 
ttik  =  ±1/2  axe  the  azimuthal  quantum  numbers  of  the 
Rb  electron  spin  and  the  Xe  nuclear  spin,  correspond¬ 
ingly.  These  are  eigenstates  of  the  zeroth-order  Hamil¬ 
tonian,  =  (hums  ±  otmsmK)\msmK)  « 

hujms\msmK)-  If  the  molecule  is  formed  in  the  states 
1 1/2, 1/2)  or  |  - 1/2,  —1/2)  there  will  be  no  spin  exchange 
since  =  0  for  these  states.  The  remain¬ 

ing  two  states,  which  we  denote  with  |1)  for  ms  =  1/2, 
mK  —  -1/2  and  |2)  for  ms  =  -1/2,  mj<  =  1/2  are  cou¬ 
pled,  and  one  can  readily  verify  that  H$  —  (1\H^\2)  = 
a/2.  From  first-order  perturbation  theory,  the  proba¬ 
bility  that  a  molecule,  formed  in  the  state  |1)  at  time 
t  =  0,  will  have  made  a  transition  to  the  state  |2)  at  time 
t  >  0  is  W  =  (a/fiu;)2sin2(a;t/2).  Previous  experimental 
work  [15,  30]  has  shown  that  XeRb  molecules  are  bro¬ 
ken  up  by  collisions  with  other  xenon  atoms  or  nitrogen 
molecules.  The  mean  lifetime  r  was  found  to  be  very 


nearly  1/r  =  ([N2]  +  3.6  [Xe])  x  5.5  x  109  s"1  amg"1.  A 
Xe  atom  is  approximately  3.6  times  as  likely  to  break  up 
a  XeRb  molecule  as  an  N2  molecule.  For  cell  T094  we 
find  r  =  31  ps.  Therefore  ru  =  51.  The  magnetic  field 
is  large  enough  that  the  probability  of  a  spin  flip  from 
state  1 1)  to  1 2)  is  independent  of  the  molecular  lifetime  r 
and  is  (W)  =  (a/tko)2/ 2  =  7.3  x  10-9.  The  factor  of  1/2 
is  from  the  mean  value  of  sin2(cjt/2). 

A  XeRb  molecule  can  form  in  simultaneous  collisions 
of  a  Xe  atom,  a  Rb  atom  and  a  third  atom  or  molecule. 
The  third  body  is  another  Xe  atom  or  a  nitrogen  molecule 
N2  in  our  cells.  The  number  density  of  molecules  will 
therefore  grow  at  the  rate 

4  [XeRb]  =C[Xe][Rb].  (A2) 

at 

where  £  =  Zr  ( [N2]  +  3.6  [Xe] ).  Three-body  rate  coef¬ 
ficients  Zr  have  been  determined  from  previous  experi¬ 
mental  work.  For  nitrogen  [15],  Zr  =  3.9  x  10“32  cm6/s 
(after  rescaling  the  original  data  by  the  ratio  of  tem¬ 
peratures  349K/453K).  For  cell  T094  this  implies  that 
(  =  6.2  x  10" 12  cm3/s.  Denoting  the  number  densities 
of  spin-up  and  spin-down  atoms  with  the  supscripts  + 
and  —  respectively,  we  see  the  the  formation  of  XeRb 
molecules  will  cause  [Xe+]  to  change  at  the  rate 

^[Xe+]  =  C([Xe_][Rb+]  -  [Xe+][Rb_])(W0.  (A3) 

We  have  neglected  the  slight  difference  in  the  forward 
and  reverse  rates,  an  effect  we  will  analyze  in  the  next 
appendix.  In  keeping  with  this  approximation  we  have 
neglected  the  thermal  polarization  of  the  Rb  atoms.  Set¬ 
ting  [Rb±]  =  [Rb]/2,  and  [Xe±]  =  [Xe](l/2  ±  (Kz))  in 
(A3)  we  find 

±{Kz)  =  -k'  [Rb](Kz),  (A4) 

where  the  three-body  rate  coefficient  for  cell  T094  is 

k  =  (  (W)  =  4.5  x  10~2Ocm3/s.  (A5) 

This  is  3900  times  smaller  than  the  experimental  re¬ 
sult,  k  =  1.75  x  10" 16  em3/s.  Similarly  small  estimates 
are  found  for  the  other  cells.  We  conclude  that  XeRb 
molecules  cannot  contribute  more  than  0.1%  to  the  rates 
measured  in  these  experiments.  The  estimated  XeRb 
contribution  is  substantially  smaller  than  the  experimen¬ 
tal  error.  The  uncertainties  in  the  estimates  of  the  molec¬ 
ular  formation  and  breakup  rates  are  not  large  enough  to 
change  this  conclusion.  Furthermore,  our  experimentally 
measured  rate  coefficients  have  no  systematic  dependence 
on  the  gas  pressure  of  the  cells,  and  are  the  same  to  within 
the  experimental  error.  This  is  to  be  expected  if  the  rates 
are  due  to  binary  collisions,  but  not  if  the  rates  include  a 
substantial  contribution  from  the  formation  and  breakup 
of  XeRb  molecules. 
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APPENDIX  B:  FORWARD  AND  REVERSE 
RATES 

Here  we  discuss  the  meaning  of  “the  binary  rate  coef¬ 
ficient”  k,  which  was  measured  in  these  experiments.  A 
spin-flip  binary  collision  where  the  electron  spin  flips  up 
while  the  nuclear  spin  flips  down  is  endothermic  by  the 
sum  of  the  electron  and  nuclear  spin  flip  energies, 

AE  =  huj  +  2iiKB.  (Bl) 

The  magnetic  moment  of  the  129Xe  nucleus  is  /xjf .  We 
denote  the  rate  coefficient  for  such  endothermic  collisions 
by  Ki2-  The  inverse  collision,  in  which  the  electron  flips 
down  while  the  nucleus  flips  up,  is  exothermic  by  A E. 
We  will  denote  exothermic  rate  coefficient  by  «2i.  De¬ 
tailed  balance  requires  that  these  two  rates  differ  by  a 
Boltzmann  factor.  They  can  therefore  be  written  as 

K21  =  KeAE/2kT,  and  k\2  =  Ke~  AE/2kT.  (B2) 

where 

K  =  y/Kl 2^21  -  (B3) 

For  binary  collisions,  the  rate  equation  analogous  to 
(A3)  is 

4[Xe+]  =  K12[Xe_][Rb+]  -  K21[Xe+][Rb_l.  (B4) 

at 

Inserting  [Rb±]  =  [Rb](l/2  ±  (Sz))  and  [Xe±]  = 
[Xe](l/2±  (Kz))  for  the  number  densities  of  spin-up  and 
spin-down  atoms  into  (B4)  we  find 

jt(Kz)  =  [Rb]fc21(l/2  —  (Kz))(l/2  +  (Sz)) 

-  [Rb]«12(l/2  +  (Kz))(  1/2  -  (Sz)).  (B5) 


We  assume  the  transition  rate  between  the  two  electron 
spin  states  of  the  Rb  atom  due  to  spin  exchange  with  the 
nuclei  of  129Xe  atoms  is  much  slower  than  the  transition 
rate  due  to  spin-rotation  collisions,  diffusion  to  the  cell 
walls,  and  other  relaxation  mechanisms  independent  of 
the  nuclear  spin  polarization  of  the  129Xe.  Then  we  can 
take  the  electron  spin  polariation  in  (B5)  to  be  time- 
independent  and  equal  to  the  thermal-equilibrium  value, 


(SZ)T  =  (1/2)  tanh(ftu;/2/cT).  (B6) 


Substituting  (B6),  (B2),  and  (Bl)  into  (B5)  we  find 


jt(Kz)  =  ^((Kz)T-(Kz)),  (B7) 


where  ( KZ)T  =  (1/2)  tanh(/x/c^/^)  »  Ri<B/2kT.  The 
solution  of  (B7)  is  Eq.  (1).  The  longitudinal  relaxation 
rate  is 


1  _  cosh^K  B/kT)  =  0  9996  [Rbj  (B8) 

Ti  cosh(fhL>/2kT)  1  J  1  J  v  ' 


The  ratio  of  the  cosh  functions  is  0.9996  for  B  —  9.4  T 
and  T  =  450  K,  the  mean  temperature  of  our  experi¬ 
ments.  So  the  measured  slope  of  1/Ti  versus  [Rb]  gives 
k,  the  mean  of  the  endothermic  and  exothermic  rate  co¬ 
efficients,  to  better  than  one  part  per  thousand. 
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